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We demonstrate the cooling of a two species ion crystal consisting of one 9 Be + and one 24 Mg + ion. 
Since the respective cooling transitions of these two species are separated by more than 30 nm, laser 
manipulation of one ion has negligible effect on the other even when the ions are not individually 
addressed. As such this is a useful system for re-initializing the motional state in an ion trap 
quantum computer without affecting the qubit information. Additionally, we have found that the 
mass difference between ions enables a novel method for detecting and subsequently eliminating the 
effects of radio frequency (RF) micro-motion. 
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I. INTRODUCTION 

A promising system for the development of a quan- 
tum computer is a collection of cold trapped ions. In 
this scheme, information is stored in the internal states 
of the ions, and logic gates are performed by coupling 
qubits through a motional degree of freedom. The orig- 
inal proposal by Cirac and Zoller P, 0] requires the sys- 
tem to be initialized in the motional ground state, with 
imperfect ground-state occupation resulting in a loss of 
gate fidelity. Other gate implementations that relax this 
condition have been proposed [3,0, HUES 111 anc ^ demon- 
strated [ij El but many of these schemes [j, El pi 13 re- 
quire the ions to be in the Lamb-Dicke limit Main- 
taining this condition in a large scale device places strin- 
gent requirements on allowable heating rates. Further- 
more, one proposed architecture for a large-scale device 
[llL requires separation and shuttling of ions between 
different trapping regions in an array of interconnected 
traps and initial experiments reported substantial heat- 
ing during the separation process |l3j . For these reasons 
it is expected that cooling between gate operations will 
be needed in a viable large scale processor. 

Ion cooling is typically achieved by laser cooling, which 
requires internal state relaxation. Therefore, direct laser 
cooling of the qubit ions is not possible without destroy- 
ing the coherence of the qubit state. Alternatively, ad- 
ditional refrigerant ions can be laser cooled directly, al- 
lowing the motional degrees of freedom to be sympathet- 
ically cooled via the Coulomb interaction [l4|. For this 
strategy to work, the cooling radiation must not couple 
to the qubit's internal state. Thus the cooling radiation 
must be sufficiently focused onto the refrigerants and/or 
be sufficiently detuned from transitions within the qubit's 
internal state manifold. In the context of quantum com- 
puting, there have been two previous implementations of 
this scheme. In one approach the refrigerant ions were of 
the same species as the qubit ions, and successful sym- 
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pathetic cooling hinged on the ability to individually ad- 
dress the trapped ions [l^. While this implementation 
was able to achieve 99% ground state cooling, individual 
addressing can be technically demanding, particularly in 
tightly confining traps, which are required to maintain 
fast gate speeds 0, El ■ In a second approach the re- 
frigerant ions were differing isotopes of the same atomic 
species Since isotope shifts of the relevant cool- 

ing transitions are typically several gigahertz and much 
larger than natural linewidth, the required degree of in- 
dividual addressing is reduced. Furthermore, the isotope 
shifts can be small enough so that optical modulators 
can provide the cooling light without the need for addi- 
tional lasers However, not all atomic species have 
isotopes that would be suitable for sympathetic cooling. 
Moreover, the amount of sympathetic cooling required 
may require cooling transitions that are much further 
detuned from any transitions in the qubit ion. For these 
reasons we have chosen an implementation in which the 
refrigerant ions are of a different atomic species, with 
the relevant dipole transitions being separated by about 
10 5 GHz (33 nm). 



The system we investigate is that of a two-ion crystal 
consisting of one 9 Be + and one 24 Mg + ion. The theoret- 
ical details of laser cooling such a crystal can be found 
in [l8L fTsj . The potential, in the small oscillation limit, 
is first expressed in normal mode coordinates. In this 
representation the system is that of six independent har- 
monic oscillators and, as such, laser cooling a single nor- 
mal mode u sing one ion is essentially the same as cooling 
a single ion [Tl| • Here we restrict our attention to the two 
modes involving axial motion only. When the masses are 
identical, these two modes are the center-of-mass (COM) 
and stretch mode, in which the ions oscillate in phase and 
180° out of phase respectively. When the masses are dif- 
ferent the two modes no longer correspond to center-of- 
mass and relative motion. However, in what follows, we 
retain the terms COM and stretch as the normal modes 
are still described by an in-phase and out-of-phase mo- 
tion analogous to the identical ion case |l9j . 
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II. SYMPATHETIC COOLING: 
EXPERIMENTAL RESULTS 

The ions are confined in a linear Paul trap in which ap- 
plied static and RF potentials provide confinement 
In the typical setup used here the axial confinement re- 
sults in COM and stretch mode frequencies of 2.05 MHz 
and 4.3 MHz respectively. A quantization axis is estab- 
lished with an applied static magnetic field of B ~ 17 G 
which is oriented at an angle of 45° with respect to the 
trap axis. For completeness we have investigated cooling 
of the two-ion crystal using either 9 Be + or 24 Mg + as the 
refrigerant ion. 

A. Beryllium Cooling 

For 9 Be + the relevant level structure is shown in 
Fig. [TJa) . Doppler cooling [2(3, 0] is achieved using the 
<t~ polarized D3 beam, with beams Dl and D2 provid- 
ing any necessary optical pumping (both of which are 
also polarized cr~). Sideband cooling is then achieved 
using beams Rl (it) and R2 (<r + /a~) to drive stimulated 
Raman transitions from \l,n) to ||,n— 1), followed by 
an optical pumping pulse provided by beams Dl and D2 
[22|. The Raman beams propagate at right angles to 
each other with R2 parallel to the quantization axis and 
the difference vector Ak = ki — k2 parallel to the trap 
axis. After 30 cooling cycles the population of the crys- 
tal's motional ground state is probed using the 9 Be + ion 
as discussed in [22, • Briefly, we drive Raman transi- 
tions on the red and blue sideband for a fixed duration 
followed by a detection pulse (beam D3), which mea- 
sures the probability of being in the \\) state. Assum- 
ing a thermal distribution over the vibrational levels, the 
ratio, r, of the red and blue sideband signal strengths 
yields a direct measure of the ground state population 
and mean vibrational quanta via P(n = 0) = 1 — r and 
n = r/(l — r) 23]. In Figs. and 01 we show cooling 
results for both the stretch and COM mode and for com- 
parison we have included results from Doppler cooling 
alone. From the sideband data we infer a ground state 
occupancy of 0.97(2) for the COM mode and 0.96(3) for 
the stretch with corresponding mean vibrational quanta 
of h — 0.03(2) and n = 0.04(3) respectively. 

B. Magnesium Cooling 

For 24 Mg + the relevant level structure is shown in 
Fig. ^b) . R2 is linearly polarized and propagates along 
the quantization axis giving equal contributions to its 
(7 + /cr~ components. This choice, together with the 7r po- 
larized Rl beam, serves to partially balance the ac Stark 
shifts induced by the Raman beams J32j. In addition, 
the intensity of the Rl beam is chosen so that it has the 
same coupling, VL (Eq. <|A2|l h to the excited states as the 
a components of the R2 beam. This choice gives the most 
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FIG. 1: Relevant energy level structures for 9 Be + and 24 Mg + 
(not to scale) (a) Level structure for 9 Be + . The 2 S\/2 levels 
are labelled by their F, tuf quantum numbers. Beams Dl, 
D2, and D3 are all polarized a~ . Raman beams Rl and R2 
are polarized n and a + ja~ respectively with a detuning A ~ 
27rx80GHz. For all beams A ~ 313 nm. (b) Level structure for 
24 Mg + . Energy levels correspond to the | J = \,mj = ±|) 
states. Dl is near resonant and polarized a~ . Raman beams 
Rl and R2 are polarized n and a + /a~ respectively, with a 
detuning A « 2tt x 750 MHz. The 17 G field defining the 
quantization axis gives a ground state Zeeman splitting of 
S w 2tt x 40 MHz. For all beams A « 280 nm. 



favorable ratio of Raman coupling to spontaneous emis- 
sion rate for this beam configuration. While a linearly 
polarized field tuned to the red of the 2 Si/ 2 2 Pi/2 
transition could provide Doppler cooling, this would re- 
quire the use of an additional laser beam. To circumvent 
this need the initial Doppler cooling step is provided us- 
ing 9 Be + as before. For sideband cooling, Raman tran- 
sitions of a fixed duration are driven by beams Rl and 
R2, which are tuned to the \gi,n) <-* \g2,n— 1) transi- 
tion and oriented similarly to the Raman beams used for 
9 Bc + . Each Raman pulse is then followed by an optical 
pumping pulse provided by the near resonant o~ polar- 
ized beam Dl. Since the 2 Si/2 «-> 2 P\/2 manifold lacks a 
closed cycling transition, we use 9 Be + as before in order 
to probe the final ground state fraction and results are 
given in Fig. From this data we infer a ground state 
occupancy of 0.84(4) for the COM mode and 0.66(3) for 
the stretch mode, with corresponding mean vibrational 
quanta of n — 0.19(6) and h — 0.52(7) respectively. 
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FIG. 2: Results for the COM mode using Be as a cooling 
ion. Red (blue) sideband data is shown on the left (right), 
(a) Doppler cooling and (b) Sideband cooling. Values of n 
are estimated from the sideband ratio - see text. 



FIG. 3: Results for the stretch mode using 9 Be + as a cooling 
ion. Red (blue) sideband data is shown on the left (right), 
(a) Doppler cooling and (b) Sideband cooling. Values of n 
are estimated from the sideband ratio - see text. 



The significant difference between the results obtained 
using 9 Be + as a refrigerant ion and those using 24 Mg + 
can be predominately attributed to spontaneous emis- 
sion from the 24 Mg + Raman cooling beams. For 9 Be + , 
the Raman beams are detuned by A m 2tt x 80 GHz and 
spontaneous emission plays a negligibly small role. How- 
ever, for 24 Mg + , all beams are derived from a single laser 
source with the Raman detuning achieved using acousto- 
optic modulators. As such our current setup limits the 
detuning of the Raman beams to A w 2ir x 750 MHz, 
and this gives rise to a significant amount of spontaneous 
emission during a Raman cooling cycle. Specifically, one 
can show that the mean number of photons scattered 
in a time corresponding to a coherent 7t pulse on the 
\g%, n = 1) «-> |<72 5 n = 0) transition is 0.55 for the COM 
mode and 2.0 for the stretch. With spontaneous emis- 
sion, population cannot be completely transferred from 
one internal ground state to the other, and the ground 
state \gi,n = 0) is no longer a dark state for the Ra- 
man cooling process. Thus, in any cooling cycle, popula- 
tion cannot be completely transferred to the vibrational 
ground state. Furthermore, after the Raman pulse, there 
will always be a finite amount of population left in the 
internal state, \g2)- This gives rise to a small amount of 
recoil heating during the final optical pumping step. 



III. SIMULATIONS 

To more firmly establish spontaneous emission as the 
principal limitation when cooling with 24 Mg + we have 
used a master equation to model the cooling process, the 
details of which are given in the appendices. In these 
simulations we have used f2 = 2tt x 30 MHz, based on 
the measured beam intensities, which yields a carrier 
Rabi frequency associated with the Raman \gi, n = 0) 
\g2,n = 0) transition of approximately 2tt x 300kHz. Ra- 
man pulse times of 2 /is and 5 [is were used, as in the ex- 
periments, for the COM and stretch mode respectively. 
The initial state was taken to be a thermal state in both 
cases with n = 4 for the COM and n = 1.7 for the stretch, 
consistent with the experimental data in Figs. 13a) and 
Ha). 

The final distributions obtained from numerical inte- 
gration of the master equation are given in Fig. [S] For 
these distributions we find ground state occupancies of 
0.80 for the COM mode and 0.77 for the stretch mode 
with corresponding mean vibrational quanta of n — 0.46 
and n = 0.31 respectively. These values of n are sig- 
nificantly different from those found in the experiment. 
However the distributions in Fig. [3] are not strictly ther- 
mal and are not well characterized by n. Recall that 
we experimentally characterize the final state based on 
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FIG. 4: Cooling results using 24 Mg + as a cooling ion. Red 
(blue) sideband data is shown on the left (right), (a) COM 
mode and (b) Stretch mode. Values of n are estimated from 
the sideband ratio - see text. 



FIG. 5: Simulation results using the master equation show- 
ing the final populations, P n , vs. the respective vibrational 
quanta, n. (a) COM mode and (b) Stretch mode. Black bars 
correspond to the simulation results and the white bars are 
the estimated thermal distributions based on sideband ratios. 



the measured sideband ratio and the assumption of a 
thermal state. Therefore, to make a better comparison 
to the experiment, we calculate the sideband ratios for 
the distributions in Fig. and we find r = 0.18 and 
r = 0.23 for the COM and stretch respectively. These 
values, together with the assumption of a thermal dis- 
tribution, give ground state occupancies of 0.82 for the 
COM mode and 0.77 for the stretch with corresponding 
mean vibrational quanta of n — 0.23 and n = 0.30 re- 
spectively. These values are in better agreement with 
experiment and, for comparison, we have included the 
estimated thermal distributions in Fig. [S] 

From the measured values of the frequencies and a nor- 
mal mode analysis we calculate Lamb-Dicke parameters 
of 0.3 for the COM mode and 0.082 for the stretch. From 
these values, together with the carrier frequency quoted 
above, one might expect pulse times of 2.8 /us and 10 /is 
for the COM and stretch respectively based on calculated 
7r times for the \g±, n = 1) <-> |<?2, n — 0) transition. The 
significant difference between these pulse times and those 
used in both the simulations and the experiments is due 
to that fact that we do not achieve ground state cooling. 
Thus the final optimum pulse need not be a 7r pulse on 
|<7i, n = 1) «-> |.92, n — 0) transition. 



IV. RF MICROMOTION COMPENSATION 

An additional technical advantage of using different 
atomic species arises from the difference in mass of the 
two ions. In a linear Paul trap [TTL |24| the transver- 
sal confinement provided by the applied RF fields can 
be described by a harmonic pseudo-potential with a fre- 
quency that is inversely proportional to the mass of the 
ion. Thus, in the presence of a transverse stray elec- 
tric field, the two ions have an equilibrium position that 
depends on the mass and their relative position vector 
no longer lies parallel to the trap axis. The net effect 
is that the modes are no longer separable into the axial 
and radial directions, and the measured frequencies devi- 
ate from those found in the absence of stray fields. Since 
this effect only occurs in the presence of stray fields, it 
provides an experimental method to detect and eliminate 
the RF-micromotion they induce. 

To determine the effectiveness of this method we com- 
pare it with a previously demonstrated method in which 
the micromotion is monitored by measuring the below 
saturation scattering rate Rq on the cycling transition 
in 9 Be + (D3 in Fig. Qfa)) when the laser is tuned near 
the carrier (wi ascr — w a tom — 0) and the scattering rate 
i?i when the laser is tuned near the first micromotion- 
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induced sideband (uvaser — w at0 m — ±^_rf) uM ■ From 
these measurements the amplitude of the micromotion 
can be found from the ratio 



R= R L= JHk-u) 



Ro J ( k ' u ) 



(1) 



where Jk are Bessel functions of the first kind, k is the 
wave vector of the laser and u is the micromotion ampli- 
tude. To compare the two techniques we determine the 
effect on the frequency spectrum for a given R. 

In the pseudo-potential approximation, the potential 
energy for a 9 Be + ion of mass m can be written 



uj\x\ + [w\ + wf) x\ + (w§ - u>l) x\ 



(2) 



where ojq is the frequency associated with the RF pseudo- 
potential and uf, u>2, and — wf are the curvatures as- 
sociated with the applied static field needed to provide 
confinement along the trap axis (xi) |33|. Since ujo and 
the static field curvatures wj? are inversely proportional 
to the mass the potential energy of the two ion system 
can be written 
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where x and y denote the positions of 9 Be + and 24 Mg + 
respectively, fi is the ratio of the 24 Mg + mass to the 9 Be + 
mass, and q is the ion's charge. In this expression we have 
included a stray electric field 
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a ^(cosy X2 + sins' X3). 



(4) 



In the limit that loq dominates the radial confinement, a 
is the off-axis displacement of 9 Be + due to the stray field. 
More generally it can be shown [25| that an electric field 
of this form leads to micromotion of amplitude 
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For our particular setup we have loq ~ 27r x 9 MHz, lo\ ~ 
2tt x 2.8 MHz, uj 2 ^ 2tt x 2 MHz, and lj 3 ~ 2?r x 3.4 MHz 
and so we neglect the small dependence of the micromo- 
tion amplitude on the applied static field. If we assume a 
probe beam is optimally aligned for micromotion detec- 
tion we have 
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FIG. 6: Plots showing effects of a transverse stray field on 
frequency spectrum for our experimental parameters, (a) 
Stretch mode shift and (b) Lamb-Dicke parameter for radial 
rocking mode in the £2 as a function of the angle 9. 



Using this expression in Eq. JQ) we can then find a 
for a given ratio R. Using this value of a in Eq. (J3J) we 
can find, for a given ft, the equilibrium position and nor- 
mal mode frequencies for the two ion system |2(j|. In 
Fig. EI a) we plot the shift in the stretch mode frequency 
as a function of 9 using R = 0.1 and the experimental 
value Qrf = 27r x 110 MHz. The maximum frequency 
shift of 90 kHz is easily detected from the difference in the 
stretch mode sideband and carrier resonant frequencies. 
The minimum frequency shift of 7 kHz on the other hand 
would indicate a reduced sensitivity to stray fields acting 
along the £2 direction. The sensitivity could be increased 
by applying a static potential with negative curvature in 
this direction. This could be accomplished by, for ex- 
ample, applying an overall static potential difference be- 
tween the RF and (RF grounded) control electrodes 0] . 
Alternatively we could utilize an additional effect that 
influences the observed Raman frequency spectrum. The 
Raman beams used to probe the mode frequencies cou- 
ple to all modes that have motion along the direction of 
the difference vector Ak = ki — k2 which, in this case, 
is aligned along the trap axis. Thus, in the presence 
of a stray field, the Raman beams can couple to other 
modes giving rise to additional features in the frequency 
spectrum. The strength of this coupling is determined 
by the Lamb-Dicke parameter, rji, and in Fig. |Sfb) we 
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plot the magnitude of 771 as a function of 9 for the radial 
rocking mode whose motion is nominally in the x 2 direc- 
tion. The maximum value of ryi ~ 0.05 gives rise to an 
easily detected feature in the frequency spectrum. Since 
this effect complements the shift of the stretch mode fre- 
quency, compensation of stray field from any direction 
can be accomplished. 

Although the sensitivity of this approach depends on 
the strength of the RF confinement, it does have ad- 
vantages over fluorescence measurements where, to null 
micromotion in the transverse direction, we require two 
non-parallel probe laser beam directions. In particular, 
monitoring the motional frequency spectrum does not 
require any additional optical access and is sensitive to 
micromotion in all transverse directions. Furthermore, 
the mode frequencies are needed for quantum logic ex- 
periments anyway and so monitoring the spectrum is a 
natural part of the experiment setup. As such, we have 
found this property to be a convenient method to detect 
and compensate the presence of stray fields. 



V. SPONTANEOUS EMISSION AND STARK 
SHIFTS 

The advantage of using a different atomic species for 
the refrigerant ion is in the small Stark shifts and off- 
resonant excitation of the qubit ion due to the presence 
of the cooling light. Although the spontaneous emission 
and Stark shifts on the qubit ion will depend on the exact 
experimental conditions, we can get an estimate of these 
effects with the following simple model. We calculate the 
probability of spontaneous emission from a qubit based 
on a superposition of the | !) and | f) states of 9 Be + , 
during a 7r pulse on the 1 01,71 = 1) <-> |<?2,7i = 0) 24 Mg + 
Raman cooling transition. (This will be the dominant 
source of spontaneous emission during one cooling cycle 
since the re-pumping pulse is assumed to be near reso- 
nant and therefore of much less intensity.) We assume 
the laser beam intensity is the same on both ions. The 
time for the cooling pulse is given by tv ~ 2ttA/ '(^il 2 ). 
From Ref. 01 j we nn d the rate of spontaneous emission 
from the 9 Be+ qubit to be R SE = (3/2) 7 ft 2 /(A*) 2 where 
A* is the frequency difference between the relevant tran- 
sitions for 9 Be + and 24 Mg + . Therefore, the probability 
of spontaneous emission from the qubit superposition is 
given by P SE = Rset* = 37r 7 A/(r/(A*) 2 ). For the COM 
mode (rj = 0.3), using the experimental values from Sec. 
Ill, we find P S e =^4x 10~ n . Better cooling with 24 Mg + 
does require us to increase the Raman beam detuning, 
which leads to a corresponding increase in the probabil- 
ity of spontaneous emission (proportional to A), but this 
should still yield a negligible probability of spontaneous 
emission. 

To estimate Stark shifts on the qubit, we first note that 
by using linearly polarized light to drive the 24 Mg + Ra- 
man coolin g tra nsition, net Stark shifts on the 9 Be + qubit 
will vanish |27| . However, during the 24 Mg + re-pumping 



pulse, we require a + light, which leads to Stark shifts on 
the qubit. From Ref. |27| . we estimate the phase shift on 
the qubit for each re-pumping pulse to be approximately 
A<fi — 3 x 10 -10 . Because the expected decoherence from 
spontaneous emission and Stark shifts are so small, we 
could not test these predictions. Nevertheless, we did 
look for decoherence due to presence of the 24 Mg+ Ra- 
man cooling beams. We performed a Ramsey resonance 
experiment on the | J.) — > | f) transition in 9 Be + with a 
precession time of 10 ms and observed no loss in contrast 
(5% accuracy) with the 24 Mg + Raman beams applied 
continuously with maximum intensity during this time. 



VI. CONCLUSION 

We have demonstrated sympathetic cooling of a two 
ion crystal consisting of one 9 Be + and one 24 Mg + . Using 
9 Be + as the cooling ion, we have been able to achieve a 
ground state occupancy for the axial modes of ~ 95% and 
better cooling results can be expected with improvements 
in the apparatus. With 24 Mg + as the cooling ion we are 
currently limited by spontaneous emission. A model has 
been developed which takes spontaneous emission into 
account and gives reasonable agreement with the exper- 
imental results. Thus, for 24 Mg + cooling, we expect 
to find a substantial improvement in the cooling results 
when we use larger detunings. For the purposes of quan- 
tum information processing one can anticipate the need 
to cool more than two ions: the three ion crystal consist- 
ing of two qubits and one refrigerant being an obvious ex- 
ample. For such cases the salient features of the cooling 
process are essentially the same as those demonstrated 
here, provided the normal mode frequencies can be spec- 
trally resolved. This spectral resolution will typically be 
achievable in the case where two adjacent qubit ions are 
cooled by a third refrigerant ion [ill IT^ . A configura- 
tion where an even number of qubit ions is cooled by a 
centrally located refrigerant ion |l2( should also be exper- 
imentally feasible. Finally, we note that these techniques 
may find use outside the realm of quantum computation, 
for instance, when applied to atomic clocks |2S[ . 

In addition to demonstrating sympathetic cooling, we 
have found that a mass difference between ions gives rise 
to a novel technique to detect and eliminate the effects 
of RF micromotion. This technique involves monitoring 
the motional frequency spectrum only and, as such, it 
is a very convenient and easily implemented method of 
micromotion compensation. 
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APPENDIX A: MODEL 

In the interaction picture the master equation has the 
usual form 



P 



-i [Hi,p] + Lp 



(Al) 



where Hi is the interaction Hamiltonian and L is the 
Liouvillian operator which accounts for dissipative pro- 
cesses. For the 24 Mg + system depicted in Fig. IHb) the 
Hamiltonian is given by 



H, 



D 



+Dt(| ei ) (g 2 \+e^ st \e 2 )(g 



i(<5 t -§<5)t 



62) (32! 
+ c.c. (A2) 



where 8l = ^112 — ^>m ~ 8 is the relative detuning of 
the Raman beams from the Raman resonance and D is 
the kick or recoil operator. In terms of the annihilation 
operator, a, the normal mode frequency, u>, and the cor- 
responding Lamb-Dickc parameter, 77, this operator may 
be written 



D 



exp 



(A3) 



where we have restricted our attention to the mode of 
interest and r\ is defined in terms of Ak for the Raman 
beams. Using the notation ajk = \gj) (e^l the Liouvillian 
has the usual form 

L P = 2 H c % (ZvjkpVjk - °\k a ikP - P a jk a Jk) (A4) 

j,k 

where Cjk is the Clebsch- Gordon coefficient connecting 
\gj) and |e^.), T is the total decay rate from the excited 
state, and p describes the density matrix after a sponta- 
neous emission event 123, 123 : 



where the upper (lower) expression applies to decay chan- 
nels involving linear (circular) polarization. 



APPENDIX B: ADIABATIC ELIMINATION 

The master equation formulated in appendix^can be 
integrated directly using a truncated basis of Fock states 
to describe the vibrational modes. However it is compu- 
tationally intensive to do so and significant simplification 
can be achieved by adiabatically eliminating the excited 
state. The procedure we adopt here closely follows that 
found in 



ic p 

m 



Let P g and P e be the projection operators defined by 
p g = 2 \9k) (9k\ and P e = J2 l e fc) ( e fc|- Tnen an Y °P er - 

k k 

ator, A, can be written as A = A ee + A eg + A ge + A gg 
where A a b — P a APf,. In this way the master Eq. (|Aip 
can be rewritten in component form giving 



Pee 
Peg 

Pgg 



-i (tl eg p ge - p eg H gi ?j - Tp e 



i [H eg Pgg — PeeH e 



1 I HgePeg PgeP-eg 



' VS )Z C )k a i,kPeeO^ 



(Bl) 
*A) /5 eg (B2) 



(B3) 



where have used the definitions Hi = e 1 H eg + c.c and 
Peg = e lAt Peg- To proceed we neglect the the terms p eg 
and Pee H eo in Eq. (|_B2f) . the validity of which is discussed 
in |3lj| . This results in an algebraic expression for p eg 
which can be rearranged to give 



dn W(k ■ e z ) exp ( ir^j=- (ae- luJt + <Je iu>t ) 



k ■ e t 



pexp I -ir]—j=- (ae luJt + a + e 



(A5) 



In this expression k is the unit vector giving the prop- 
agation direction of the scattered photon relative to the 
quantization axis denoted e 2 , St is a unit vector along 
the trap axis, W(k ■ e z ) is the angular distribution of the 
emission, and the integration is carried over the full solid 
angle. For the dipole decay discussed here we have 



(A6) 




Peg = 



(r/2) + lA 



Heg Pgg- 



(B4) 



Substituting this expression into Eq. I|B1|I and neglecting 
p ee then yields the result 



Pee = 



(r/2) 2 + a- 



-H eg p gg H ge . 



(B5) 



Finally, Eqs. (|B4|> and (|BfS|) can be substituted into 
Eq. (|B3|I to give a closed equation of motion for p gg . 
Dropping the subscripts on p gg we find the equation 



8 



-A 



{T/2f 



2 H ge H eg , p 



(r/2) 



(r/2) 2 + A 2HgeHeg,p 
i 



(r/2) 



(r/2) 2 



c jk 2uj k H eg 



(B6) 



for the effective ground state density matrix, where { , } 
is used to denote an anti-commutator. If one makes the 
identifications 



H eff z 

r' = 



:H„ K H K 



(r/2) 2 + a 2 9e e9 ' 

T 



(r/2) 2 + A 2 ' 



and 



jkfl eg 

then Eq. I|ij6l) can be recast into the usual form given in 
Eqs. IjAlfl and l|A4(l . Numerical simulation of the reduced 
master equation can be further simplified by expanding 
out the terms in (|B6|) and neglecting the off-resonant 
terms associated with the Raman pair er + and 7r. 

The resulting system gives excellent agreement with 
the full system given in Eq. (|A1|) . However the proce- 



dure outlined here is not rigorously correct due the resid- 
ual time-dependence in H eg associated with the Zeeman 
splitting of the excited and ground state manifolds. In 
effect, the procedure here neglects the small differential 
Stark shift induced by this splitting as well as the small 
change in the various spontaneous emission rates. These 
have been accounted for in a more elaborate treatment 
in which the projection operators are split allowing the 
individual level shifts to be accounted for. The resulting 
algebra is more involved but the procedure is precisely 
the same. For brevity and clarity we have omitted the 
details of this more elaborate treatment. Einally, we note 
that the optical pumping process can be treated in a sim- 
ilar manner with Eq. (|A2|) replaced with 



(B7) 



Hi = -D^\e x ) (g 2 \ + c.c. 
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